Supercurvature Modes from Preheating in an Open Universe 
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Preheating following the inflationary phase in models of open inflation is considered. The most 
significant difference from a preheating scenario in flat space is that supercurvature modes can be 
populated and amplified in the open universe. In certain models, such modes can dominate the usual 
resonant particle production, thus altering the ensuing thermal history within an open universe. The 
accuracy to which the masses and couplings need to be tuned to produce such supercurvature modes, 
however, makes any large deviation from the flat space cases unlikely. 
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I. INTRODUCTION 

There has been much activity lately in two exciting ar- 
eas of inflationary cosmology: the construction of "open 
inflation" models, and the development of a new theory 
of post-inflationary reheating. To date, these two de- 
velopments have been treated separately; in this paper, 
we explore some qualitatively new behavior for the re- 
heating epoch in an open universe. The open inflation 
models exploit the fact, first noted by Coleman and 
De Luccia and independently by Gott and Statler [Q, 
that the metric on the interior of a bubble nucleated 
in a sea of de-Sitter-like false vacuum would be that 
of an open Friedmann-Robertson- Walker spacetime. By 
arranging for an epoch of "slow-roll" inflation to occur 
within the nucleated bubble, and arranging for the num- 
ber of e-folds of expansion during this second inflationary 
phase not to exceed N ~ 65, it is possible to produce an 
observable universe whose mass density today falls be- 
low the critical mass density required for a flat universe: 
£! Q ~ 0.2 — 0.3, rather than the standard inflationary pre- 
diction, fl — 1.0. 

Independently, several groups have re-examined re- 
cently the process of particle production at the close of an 
inflationary phase, and have found resonances and expo- 
nential instabilities in the decay of an inflaton field into 
bosons, which had been overlooked in the original ac- 
counts of reheating. [TT|| According to the new theory 
of reheating, a post-inflationary second-order phase tran- 
sition should be viewed as having three distinct stages, 
rather than the older two-stage picture: the highly effi- 
cient, often resonant, decay of the oscillating inflaton field 
into boson decay products (which might include inflaton 
bosons themselves, due to the inflaton's self- interaction) , 
followed by the decay of these decay products (along the 
methods originally developed in frL2|| ), followed finally by 
the thermalization of these particles and the establish- 
ment of the radiation-dominated era. Kofman, Lindc, 
and Starobinsky term this first, resonant stage the "pre- 
heating" stage. H 

As demonstrated in this paper, under certain condi- 
tions the resonant decay of the inflaton in an open uni- 



verse can populate "supercurvature" modes, which have 
no analogue in the flat space case, and no particle-like 
interpretation. In some open inflation models the am- 
plification of such supercurvature modes can even domi- 
nate the usual resonant creation of particles at the time 
of reheating. Thus, models of open inflation can display 
dramatically different thermal histories following the end 
of inflation from inflationary scenarios in a spatially flat 
universe. 

A Friedmann-Robertson- Walker universe with nonva- 
nishing spatial curvature has a physical curvature scale 
given by a(t)/\K\, where a(t) is the cosmic scale factor, 
and K = — 1 or +1 for an open or closed universe, re- 
spectively. With K = — 1, the comoving curvature scale 
is thus simply +1. From the Friedmann equation it can 
further be shown that when K = — 1 the physical cur- 
vature scale will always be larger than the Hubble dis- 
tance, for any cosmological epoch. Eigenfunctions 
of the Laplacian operator can be found for this back- 
ground spacetime with eigenvalue — {k/a) 2 , where (k/a) 
is the inverse of a physical length: k/a = 2n/\ p i l , with 
< k 2 < oo. Modes with comoving wavenumber k 2 > 1 
thus vary on scales shorter than the comoving curvature 
scale, and are labelled "subcurvature" modes, whereas 
modes with < k 2 < 1 correspond to "supercurvature" 
modes. (See, e.g., |L3|].) (Note that tensor modes can 
have different values from scalar modes, such as k 2 = —3 
for gravitational waves. @ 0) In a flat universe, the 
comoving curvature scale runs off to infinity, so there are 
never supercurvature modes. It is this difference from 
the flat space case that we examine here, for the epoch 
of preheating: even though fl — * 1 at the end of the sec- 
ond round of inflation, there will still exist a comoving 
curvature length scale at +1, the existence of which may 
change the usual physics at the time of preheating. 

We begin the study following the nucleation of a sin- 
gle bubble in the midst of a sea of de Sitter space false 
vacuum. Inside the nucleated bubble, there is a nonva- 
nishing scalar field potential and a single relevant scalar 
field to drive a second round of inflation within the bub- 
ble; both the single-field models of open inflation [|l 
and Linde's two-field "hybrid" open inflation model 
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fit such a description. Inside the bubble, the Lagrangian 
density can thus be written: 
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R 



^g" u d^pd v p 



V(<p) 



(1) 



and the background spacetime takes the form Q 
ds 2 = -dt 2 + a 2 (t) [d X 2 + sinh 2 X {dO 2 + sin 2 9d<j) 2 ) 



(2) 



If we make the familiar decomposition, 

<p(t, x, 0, <f>) = <p (t) + S4>(t, X , 9, (f>), (3) 



then the equation for the quantum fluctuations 6<p be- 
comes 



□ (<ty) - 



f d 2 v 
\d<p 2 



(Stf,) = 0. 



(4) 



Switching to conformal time, drj = a 1 dt 1 and us- 
ing the expression for the d'Alembertian, Dip = 
(-gy^dJi-gy^g^dvip], then inside the bubble 
equation (Eh takes the form: 



d 2 2 a' d 
drj 2 a drj 



d 2 v 

dp 2 



(5) 



where primes denote derivatives with respect to con- 
formal time, rj, and L 2 is the comoving Laplacian op- 
erator. Eigenfunctions of L 2 in the metric of equa- 
tion (2) have been studied by many people (see, e.g. 
Jl6],[l7[). Denoting these harmonics as Yfci m (x? 9, (/)), with 
]~> 2 Ykim = —k 2 Yki m , then we may parametrize the fluctu- 
ations in terms of mode functions as: 



*>k(>7> X, 
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a(r?) 



i>(r))Ykim(x> °, 4>). 



(6) 



(We will not need the explicit form of the Yfeim's here, 
which can be written in closed form in terms of asso- 
ciated Legendre polynomials. fl6}| ) Finally, if we write 
(d 2 V/dip 2 ) iPo as M 2 ((p ), then the equation for the infla- 
ton fluctuations becomes: 



Vi' ipu + k 2 ^ k + a 2 M 2 ip k = 0. 

a 



(7) 



(We have neglected the metric perturbations which cou- 
ple to tpk , which should be a good approximation during 
the reheating epoch; for more on such metric perturba- 
tions, sec |18| .) In this paper, we will be concerned with 
the behavior of solutions of equation (^) after the end of 
inflation, as the inflaton field p oscillates near the min- 
imum of its potential. 

In Section II, we briefly review the conditions under 
which the oscillating inflaton may decay resonantly into 



bosons. In an open universe, with the possibility of pop- 
ulating supercurvature modes during the inflaton's para- 
metric resonance, we must pay closer attention than is 
usually done to the conditions under which the inflaton's 
decay may be resonant: the harmonic relations between 
the frequency of the inflaton's oscillation and the decay 
products' momenta are specified here for all resonance 
bands within the narrow resonance regime. In Section 
III, these resonance conditions are studied for the sim- 
plest model of an inflaton decaying into inflaton bosons, 
due to a p A self-coupling. As will be shown, in this 
case no supercurvature modes may be populated dur- 
ing preheating, and there should be no deviation from 
the usual flat space case. The model is generalized in 
Section IV to the case of an inflaton decaying into a dis- 
tinct species of boson, via Yukawa and quartic couplings, 
and in these cases supercurvature modes can be popu- 
lated, if the masses and couplings of the model satisfy 
certain relationships to very high accuracy. None of the 
resonant behavior will affect inflaton decay into fermions, 
because of Fermi-Dirac statistics, and we will ignore such 
couplings here. The results are discussed in Section V. 



II. PREHEATING AND RESONANT DECAYS 

Near the minimum of its potential, the inflaton will 
begin to oscillate, and this can set up a parametric res- 
onance in the inflaton's decay to bosons, j^j Assum- 
ing that the frequency of the inflaton's oscillations is 
greater than the Hubble parameter at this time, we 
may neglect Hubble expansion during the course of the 
resonant decays. Then the inflaton field oscillates as 
<Po(v) — Wo cos(w sc?7), where ~pZ is slowly decreasing on 
the oscillation time-scale = (a e m) _1 . Here a e is a 
constant giving the value of the cosmic scale factor at 
the end of inflation. Note that (a e i/ e ) _1 < 1, as will 
be true for any cosmological epoch in an open universe. 
Hence, our assumption about the frequency of the infla- 
ton's oscillation translates to (a e m) _1 (a e i? e ) _1 < 1, 
or a e m ^> 1. The A4 2 (p ) term then contains two types 
of terms, oscillating and non-oscillating, and may be writ- 
ten: 



M 2 {tp ) = A + B cos(gw osc ry), 



(8) 



where A, B, and q are determined by a given model's 
particular couplings. The equation of motion for the in- 
flaton fluctuations, equation ([?]), then becomes: 



i>k + u l [1 + 9 cos(quj osc ri)} ip k ~ 0, 



to 
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k 2 + a 2 e A , g = a 2 e B/uj 2 k . 



(9) 



This is now in the form of the Mathieu equation, solu- 
tions of which reveal exponential instabilities within var- 
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ious resonance bands .0 

The narrow resonance regime of equation (|^), with 
g < 1, can be studied perturbatively. The analytic anal- 
ysis in dJlil demonstrates that in this case the solutions 
are exponentially unstable in narrow resonance bands 
nquj asc = 2idk ± \ with |e n | < ojk, and n a positive 
integer. Each resonance band with n > 1 corresponds to 
keeping terms of 0(g n ) in the solution to equation (||), 
and the width of each successive resonance band in the 
narrow resonance regime shrinks as fllEj] : 



n 2n - 3 g n u> k 
2 3(n-l) 1)! 



= Kg n u k . 



(10) 



Near the center of each resonance band, solutions take 
the form 



rk', 



to) 



exp(±4 7?) 



with real and positive. Modes with a given k 2 
out of the resonance band due to the slow decrease of 
Tp^, and hence of g, due to scatterings and back-reaction 
effects. Given such a time-dependence of g, modes will 
only become resonantly amplified if they initially satisfy 
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nqoj 



2w k 



1, 



(12) 



modes initially satisfying nqta osc — 2ui k — f l e n|j on the 
bottom edge of the resonance band, will not remain 
within the band long enough to become amplified. We 
will thus study equation (|l^) for various models to deter- 
mine which areas of parameter space allow for the pro- 
duction and amplification of supercurvature modes. The 
subcurvature modes should behave as in the flat space 
case, and hence we will not treat them further here. 

As demonstrated explicitly in B, including a nonzero 
Hubble expansion during the period over which the in- 
flaton oscillates does not change the qualitative form of 
these solutions, but merely narrows the width of the res- 
onance bands. Similarly, including back-reaction and 
rescattering effects should further decrease the width of 
each resonance band. |||,^|J|JT^,^,^l| We may accom- 
modate these effects by inserting a phenomenological 
parameter a(j]), with < a(ri) < 1, into the band 
width, writing the width of each resonance band as 
l^nl = b n a n (r])g n ujk, though we will set a = 1 in the 
following. 

The broad resonance regime of equation (^) (gj§> 1), 
identified in || , has been studied analytically in |]10| , and 



*As emphasized in JqJlQ], the exponential instabilities arise 
for solutions of any differential equation with periodic coeffi- 
cients, and are more general than the particular form of equa- 
tion (|. 



also leads to solutions of the form ipt 



with 



» 



real and positive (and a k 



(n) 



tx exp(±cr^ ) ry), 
> 4 n) )- How- 



ever, such a broad resonance regime is ruled out for most 
quartic couplings fici[| ; and even for models with Yukawa 
couplings between the inflaton and a distinct species of 
boson, the coupling strength (and hence g) may not 
be made arbitrarily large, because such large couplings 
would interrupt the usual inflationary dynamics prior to 
the preheating epoch. Furthermore, the particles pro- 
duced during a broad resonance have much higher mo- 
menta than particles produced in the narrow resonance 
regime, and hence no supercurvature modes are expected 
to be populated from a broad resonance. For these rea- 
sons, we will restrict attention here to the narrow reso- 
nance regime for each of the models studied below. 



III. INFLATON DECAY INTO INFLATON 
BOSONS 



In this section, we consider a model in which the oscil- 
lating inflaton field decays into inflaton bosons, due to a 
<p 4 self-coupling. Near the minimum of its potential, any 
model of open inflation should assume a general chaotic 
inflation form [E2|: 



(13) 



In this case, the parameters in equation ([)]) become 
uj 2 k = k 2 + aim 2 + \a 2 e \Tpf , g = ^a 2 e \Tpf / 'ui 2 k , (14) 



or, combining terms, 

2 _ (k 2 + a\m 2 ) 



0-9) 



(15) 



Note from the form of g in equation (14) that g < 1 



always in this model: in the case of inflaton decay into 
inflaton bosons, the only possible resonances lie in the 
narrow resonance regime, and hence the explicit expres- 
sions in section II may be used. Finally, for the potential 
of equation (|l3|), the oscillating term in equation ([)]) os- 
cillates at frequency 2lo osc = 2a e m, or q = 2, after use 
has been made of a double-angle formula for cos 2 (tu osc r)) . 

Given that a e m 1, we need only consider wavenum- 
bers in the range < k <g; a e m in our search for any 
supercurvature modes which may be resonantly ampli- 
fied during preheating. Defining 



(16) 



a e m, 



then equation ( |l2] ) becomes: 



2 [n^/l 
1 



0(£ 6 ) 



1 - - A K9 n 



1 + f n g n 



(17) 
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Keeping only the lowest-order terms in £, this gives 



2alm 2 



bn9 n 



4 + b n g n 



(18) 



When n = 1, the quantity on the righthand side of 



equation (18) is negative, revealing that neither super- 
curvature nor subcurvature modes will be amplified in 
the lowest resonance band when k < a e m. For n > 1, 
the approximation £ <§; 1 is violated, and thus equa- 
tion ( Jig ) cannot be used to study resonant subcurvature 
modes. Instead, it can easily be demonstrated that when 
£ > 1 (and hence for subcurvature modes only), resonant 
modes correspond to 



ny/l-g 
2 7 



i± Wi 



2 7 2 



n 2 (l-<?) 



7=1 + -b n g n > 1, 



0(fc^ 2 ), 
(19) 



again revealing that no consistent solution (with fc res real 
and positive) may be found for n = 1. Thus, for the po- 
tential of equation (|l3|), only subcurvature modes may 
be amplified during preheating, and only in the narrow 
resonance regime with n > 2. 

For this particular potential, a limit on g may be set 
by the independent constraints on m and A from ob- 
served microwave background anisotropies: A < 10 - , 
and m/Mpi < 1CT 6 . |§ An upper bound may be placed 
on TpZ by equating this with the value of tp once inflation 
ends, that is, once the slow-roll conditions are violated, at 
d 2 V/dip 2 ~ 24nM~ l 2 V. This gives ~ < 0.16M p i (using 
the constraints on A and m), or \TpZ 2 /m 2 < 2.6 x 10~ 4 , 
guaranteeing that g < 4 x 10~ 4 . With so small a g, res- 
onance bands with n > 2 are unlikely to produce large 
occupation numbers of inflaton bosons, and the resonant 
decay of the inflaton may be rather inefficient in trans- 
ferring the energy density of the decaying inflaton field 
into boson decay products. Such tight constraints on 
the potential's parameters are absent in the case of infla- 
ton decay into a distinct species of boson, assuming that 
the second boson field is unimportant during the infla- 
tionary phase (when the primordial density perturbation 
spectrum is determined). It is to these models that we 
now turn. 



IV. INFLATON DECAY INTO A DISTINCT 
BOSON SPECIES 

In the case of inflaton decays into a distinct species of 
boson, x, we may study a potential of the form: 



In supersymmetric models, these couplings will satisfy 
H = 2\ x m. m Writing X = a" 1 ^)/^)^™^), the 



-> d 2 V/d x 2 = m\ + 2 W + 2\ 2 x ip 2 . 
(a e m) _1 <C (agi/)" 1 , the ^-modes 



equation for f(rf) will assume the same form as equa- 
tion (@), with M 
Again taking uj~^ c 
will obey: 



j'l + uj 2 x [1 + hi cos(2cj osc ?7) + h 2 cos(cj OS c?7)] fk - 0, 
lu 2 x = k 2 + a 2 e m 2 x + a 2 e \ 2 x TpZ 2 , 

u — 2 \ 2 2 / 2 t _ o 2 ; 2 

hi = a e \ tp /u , h 2 = 2a e Li(p /uj 



(21) 



Note from the form of hi that hi < 1 always, so only the 
narrow resonance regime is viable for the quartic cou- 
pling. If we first consider h 2 — > 0, and if the mass dif- 
ference between the inflaton and the x boson is small, 
then modes with k <C a e m x will become exponentially 
amplified in bands centered on: 



2a\m 2 x 



An(m/m x )^l — hi 



b n h\ 



(22) 



The next order term is of 0(k A a e 2 m x 2 ), which will be 
completely negligible if < k 2 < 1 and a e m x ^> 1. 
Though the form of equation ^22\ ) for k 2 cs looks quite sim- 
ilar to the corresponding expression for ip — > tp decays, 
equation jlq), there are important differences. First, 
there are now two free parameters, (m/m x ) and hi, 
rather than only the single free parameter g in equa- 
tion (p"8|). Second, the coupling constant A x (and hence 
hi) is not limited directly by today's observed microwave 
background spectrum, because it is assumed that the \- 
field plays no role during inflation, so hi need not be 
constrained to hi ~ O(10~ 4 ) as g is. Finally, the ra- 
tio (m/my) may be either greater or less than unity: as 
noted in , it is possible for the inflaton to decay into a 
particle heavier than itself during the resonant preheating 
stage (a decay which is kinematically forbidden accord- 
ing to non-resonant Born theory). Thus, supercurvature 
modes will be produced and resonantly amplified in a 
resonance band n* which satisfies 



< 



Sn*(m/m x )v / l — hi — 8 — 2b 1lt h"" 



< (a e m x Y 



(23) 



If a e ~ e 65 ~ 10 28 , then this is a very tiny window: the 
resonant production of supercurvature modes requires 
that the masses and couplings of the potential obey cer- 
tain interrelationships to exponential accuracy. Still, by 
exploiting the freedom in (m/m x ) and A x , it is possible 
to arrange for such conditions to be met (and the nearly- 
exact cancellation of combinations of parameters is not in 
itself foreign to particle physics). Furthermore, if n* = 1 
and hi -C 1, then all higher (subcurvature) modes could 
fall outside of resonance bands, leaving only the super- 
curvature modes to be produced and amplified at the 
preheating stage. Note that if m x <C m, then equation 
( p2| ) is inappropriate, and only subcurvature modes will 
be amplified at preheating. Other resonant subcurvature 
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modes, with k > a e m x , may be found analogously to 



equation (|19|), and will not be pursued here. 

Similarly, if hi — > and the Yukawa coupling in equa- 
tion ( ^p|) is active (hi ^ 0), then supercurvature modes 
will be produced resonantly in a resonance band if 



< 



4n*(m/m x ) — 8 — 26 n ,^2 



< (a e m x Y 



(24) 



the narrow resonance regime, up to terms of 



0{k 4 a e 



TO, 



"). Again, the pair of parameters (m/m x 



and \i may be tuned to satisfy this condition, and if 
?i* = 1 and h<i -C 1, then the production and amplifica- 
tion of supercurvature modes will dominate the preheat- 
ing phase. As for the h\ case, if m x <C m or if k > a e m x , 
then only subcurvature modes will be produced at pre- 
heating. 

Care must thus be taken to identify the (model- 
dependent) dominant decay channel of the inflaton, es- 
pecially if the two couplings in equation (ECN) are in com- 



petition with the inflaton self-coupling in equation (13). 
The conditions under which the supercurvature modes 
may become amplified and dominate the production of 
subcurvature modes during the preheating phase depends 
sensitively upon the couplings and mass scales of the de- 
caying inflaton and its decay-product bosons. Note that 
although fc^ es —>■ for these supercurvature modes, there 
is no infrared catastrophe, since every resonance band in 
the narrow resonance regime has a finite band width. 



V. DISCUSSION 

What do these exponentially amplified supercurvature 
modes represent? In the usual flat space preheating sce- 
nario, narrow resonance modes simply correspond to the 
production of low-momentum (small k) particles, whose 
occupation numbers may be estimated with the aid of a 
Bogolyubov transformation. Yet in the open universe, 
such small-/c modes can stretch beyond the curvature 
scale (and hence beyond the horizon), and no particle 
interpretation may be given to them. (This is a specific 
example of how tenuous the notion o f "p article" can be- 
come in general curved spacetimes. |23f|)P1 Instead, they 
will contribute to (very) long-wavelength curvature per- 
turbations, akin to supercurvature modes produced dur- 
ing inflation ]T^-[T^| , though their evolution will be com- 
plicated by the highly non-adiabatic nature of the re- 
heating epoch. The detailed evolution of these modes 



may be tracked using the methods of |18|, though it 
is unlikely that they will contribute greatly to any ob- 
served microwave background anisotropies today, given 
their very long correlation lengths [Q; this is a subject 
of further study. 

Rather, their most dramatic role could be in chang- 
ing the thermal history immediately following reheating: 
the new reheating scenario depends upon the explosive 
production of particles, far from thermal equilibrium; 
these decay products then thermalize via interactions 
on a (potentially) quite different time-scale. When the 
occupation numbers of out-of-equilibrium decay-product 
particles becomes quite large, this opens up the possi- 
bility for such microphysical processes as non-thermal 
phase transitions |24j| , preheating- induced supersymme- 
try breaking |2f|, and GUT-scale baryogenesis p6| . Yet 
the foregoing analysis reveals that in an open universe, 
the (resonant) preheating production of particles may be 
either totally absent or strongly diluted, due to some por- 
tion of the inflaton's energy density being 'siphoned off' 
into supercurvature modes. For models in which the res- 
onant production of subcurvature modes is subdominant, 
the inflaton will still decay into particles and populate the 
universe following the end of the resonance regime, ac- 
cording to the processes of the older theory of reheating. 
If this resonant production of supercurvature modes and 
non-resonant production of particles becomes the domi- 
nant effect, then these other microphysical processes in 
the post-preheating epoch may need to be reconsidered. 

Given the exponentially tiny supercurvature resonance 
regimes in equations (^2|) and (|24|), however, it is rather 
unlikely (though still possible) that the supercurvature 
modes from the resonant cp — > \ decays would dominate 
the preheating epoch. And such modes are completely 
absent from a strictly <p — ► (p preheating epoch. Thus, 
large deviations from preheating scenarios in flat space 
are unlikely to appear for models of open inflation. 

Finally, as noted in the conclusion in ||, for models 
with non-minimally coupled scalar fields which recover 
the Einstein- Hilbert gravitational action near the end of 
inflation (e.g., models of the form studied in (27]]), the 
analysis of reheating carries over unchanged. Thus, the 
recent open inflation model based on induced-gravity the- 
ory [^8| should display the same resonance features and 
supercurvature modes as analyzed above. 
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